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Abstract 

In this paper we consider the Chandrasekhar equations arising in optimal 
control problems for linear distributed parameter systems. The equations are 
derived via approximation theory. This approach is used to obtain existence, 
uniqueness and strong differentiability of the solutions and provides the 
basis for a convergent computation scheme for approximating feedback gain 
operators. A numerical example is presented to illustrate these ideas. 
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Introduction 


It has been noted by a number of authors [4], [7], [10], [15] that 

Chandrasekhar type algorithms can significantly reduce the computations 

necessary to calculate optimal feedback gains for linear quadratic control 
problems when the number of inputs and outputs is small, relative to the 
dimension of the state space* These algorithms were initally developed for 
finite dimensional linear time invariant systems [7], [8], [15], [16] and 

later extended to various infinite dimensional systems [4], [10], [19], 

[20]. It has been observed [4], [10], [18], [20] that the Chandrasekhar 

algorithm when applied to certain distributed parameter systems affords a 

significant computational reduction, often even greater than in the finite 
dimensional case. Therefore, it is worthwhile to consider the numerical 
aspects connected with the solution of the Chandrasekhar equations in infinite 
dimensional spaces . 

Before one can develop an approximation theory for these equations, it is 
necessary to first examine certain basic questions such as existence, 
uniqueness and regularity of solutions. The question of existence has been 
addressed by Casti and Ljung [10] for general time invariant systems 
(including certain boundary control problems) and extended to time varying 
systems by Baras and Lainiotis [4]. Sorine [18], [19], [20] developed 

existence, uniqueness and differentiability results for a special class of 
parabolic systems and noted in Reference [20] that uniqueness is much more 
difficult to obtain than existence. All of the papers cited above utilize the 
variational framework of J. L. Lions and formulated the Chandrasekhar 
equations in differential form. Distributional derivatives were used to 
define the equations and solutions which often can complicate convergence 
analysis of numerical schemes. 
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Xn the present paper we present an approach to the Chandrasekhar 
equations that is based on approximation theory. We restrict our attention to 
time invariant control systems with bounded input and bounded output 
operators. As a result we obtain existence, uniqueness and smoothness of 

solutions to integral versions of Chandrasekhar equations for distributed 
parameter control systems that include delay and hyperbolic systems not 
covered by Sorine's results [20]. Moreover, sufficient conditions for 
convergence of general approximation schemes is estabilished. 

2. The Chandrasekhar Equations 

Let H ,U and A be Hilbert spaces. We denote by L(X»/) the Banach 
space of bounded linear operators between the Hilbert spaces X an d V , 

endowed with the uniform operator topology. Throughout this paper we assume 
that A is the generator of a Cg-semigroup S(t) on H , B eL (U,H), 

V e L(H,A) and ReL (U ,U) is self-adjoint and satisfies JR# ^ m > 0. Let 

Q e L(H,H) be defined by Q = V*V where V* is the adjoint of the 

operator V. 

The linear time invariant quadratic optimal control problem is to choose 
u(») g L,(0,T;U) that minimizes the cost functional 

T 

J - J ( By ( s') B ? + <Ru(s) ,u(s)>)ds (1) 

0 

subject to the constraint that 


z(t) = Az(t) + 8u(t) 


( 2 ) 



( 3 ) 


with output 


z(0) = z Q 
y(t) = Vz(t). 


(4) 


Solutions to (2) , (3) are always defined to be mild solutions given by 

t 

z(t) = S(t)z_ + / S(t - s)8u(s)ds (5) 

U 0 

Under the assumptions stated above, it is known (see References [11], [13]) 

that there exists a unique u(») e Z_ 2 (0 ,T; (J) that minimizes J. Moreover, 
the optimal control is given by 

u(t) = -R -1 8*n(t)7(t) (6) 


where n(t) is the unique solution to the operator evolution equation 

T 

n(t) = / s*(n - t) [q - n(n)8R -1 8*n(n)]s(n - t)dn (7) 

t 

and ¥(•) is the optimal trajectory generated by TT( • ) (see References [13], 
[14] for details). It is shown in References [11], [13], [14] that the 

operator [ A - BR - ^B*II(t) ] generates an evolution operator U(t,s) 

on H and 

T 

Ji(t) = / u*(n,t)[Q + n(n)8R" 1 8*n(n)]u(n,t)dn 

t 


( 8 ) 



-4- 


A formal differentiation of equation (7) (or (8)) yields a Riccati operator 
differential equation. However, the meaning of such a differentiated form of 
the equation must be precisely defined. In the papers by Casti and Ljung [10] 
and Baras and Lainiotis [4] distributional derivatives are used to obtain 
strong solutions. Sorine [20] used a similar approach. However, Sorine was 
also able to obtain regularity of strong solutions for the special case that 
A generates an analytic semigroup. This is an important matter since the 
derivation of the Chandrasekhar equations presented in References [10] and 
[20] make heavy use of the differentiability of II( t ) . We shall avoid many of 
these difficulties by concentrating on an integral version of the 
Chandrasekhar equations. The special form of the equations can be exploited 
to obtain existence, uniqueness and some regularity properties of the gain 
operators. 

Let K(t) = R-lB*II(t) where H(t) satisfies (7) (or (8)). The goal is 
to derive a set of equations that allows one to solve for K(t) directly 
without first solving for H(t). We state here without proof the following 
theorem. 

Theorem 1 . There exists K(t) t L(H ,U) and L(t) e L(H, A) such that 
for each z e H 

T 

K(t)z = / R *B L (n)L(n)zdn (9) 

t 

T 

L(t )z = VS(T-t)z - / L(n)BK(n)S(n - t)zdn 

t 


( 10 ) 
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for 0 t T* Moreover, K(t) and L(t) = VU(T,t) are the unique strongly 
continuous solutions to (9), (10). 

The representation (9), (10) can be exploited to establish the following 

result. 


Theorem 2. If z z H, then K(t)z is continuously differentiable. 
Moreover, if z ep(A), then L(t)z is differentiable and 


^ K(t)z = -R- 1 B*L*(t)L(t)z (11) 

^ L(t)z = L(t)[A +8K(t)]z (12) 


for 0 < t <T. 

Note that (9) can also be used to establish the differentiability of the 
Riccati operator. Under the additional assumption that z e V(A) it can be 
shown that II(t)z is a continuously differentiable solution to the 
differentiated form of the Riccati equation. Although detailed proofs of 
these theorems will appear elsewhere, it is worthwhile to outline the approach 
in order to make a few observations about computational algorithms for 
approximating K(t) • 

Definition 1 . A strong approximating sequence for the control problem 

e N N N 

defined by equations (l)-(4) is a sequence (A ,V ,R J such that 
A N ,B N ,R N and Q N = [v N ] V N are bounded operators satisfying RR N H >_ m > 0, 
Q N > 0, 8 N - 8 strongly, R N + R strongly, Q N + Q strongly and the 

operators A w generate Cp-semigroups S (t) satisfying 


S N (t)z S(t)z and [s N (t) ]*z + S*(t)z 


for all z c H, uniformly for t e [0,T], 

Observe that the assumptions on (A N ,8 N ,V N ,R N ) in the definition of a 

strong approximation sequence are precisely those conditons used by Gibson 

[13], [14] (except for the boundness assumptions on A N ) to establish strong 

convergence of solutions to approximting Riccati equations to the solution 

H(t) of (7). For the case of a finite dimensional control space U, this 

strong convergence implies uniform convergence of the gain operators (see 

Theorem 6.2 in Reference [14]). However, in order to take full advantage of 

the factorization of Q = V V additional assumptions on the convergence of 
N 

V will be needed. Therefore, we state an additional condition. 

Hypothesis 1 . The operators converge strongly to V. 

Note that strong convergence of Q N to Q does not imply strong 

convergence of V N to V and, conversely. Hypothesis 1 does not imply strong 
N 

convergence of Q to Q. In most examples we have considered. Hypothesis 1 
is easily established. 

Consider the approximating system of Chandrasekhar equations 

K N (t) = / [R N ] _1 [b N ]* [L N (n)f L N (q)dn (13) 

t 

L N (t) = v^s N (T - t) - / L N (n)8 N K N (n)s N (n - t)dn 

t 


( 14 ) 
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where ( A N > B N >V N , R N ) is a strong approximating sequence. Since A N is 
bounded it is straightforward to show that (13), (14) has a unique strongly 
continuous solution. We state the following lemma. 

Lemma 1. If there exists a strong approximating sequence (A 1 *, 
such that Hypothesis 1 is satisfied, then the Chandrasekhar equations (9), 

A A 

(10) have a unique strongly continuous solution K(t), L(t) and 
K N (t) * K(t) strongly, L N (t) + L(t) strongly. 

This result can be established by a slight extension of Gibson^s results 

in Reference [13]. Lemma 1 not only provides existence and uniqueness, it 

provides sufficient conditions for the convergence of numerical approximation 

schemes. The only remaining "hole” in proof of Theorem 1 is to establish that 

there always exists strong approximating sequences that satisfy Hypothesis 1. 

Defining to be the Yosida approximation A^ = NA(NI - A) * (for N 

N N 

sufficiently large so that N e p(A)) and letting V = V, 8 =8, 

N 

R = R, the resulting sequence satisfies all the conditions of Lemma 1. It is 
important to note that Yosida approximates provide a tool for establishing 
Theorem 1 and lead to a numerical scheme that is convegent. However as a 
practical matter, Yosida approximates do not generally lead to efficient 
numerical algorithms. Therefore, in practice it is worthwhile to develop 
other approximating schemes that satisfy the sufficient conditions in Lemma 1. 
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3. A Numerical Example 

We illustrate the power of the Chandrasekhar algorithm by applying the 
averaging approximation scheme (see References [2], [3], [14]) to a simple 

delay-differential equation model of the two-dimensional airfoil shown in 
Figure 1. We note that the averaging scheme has been shown to satisfy the 
sufficient conditions for convergence (see References [2], [14]). In order to 
make efficient use of the structure of the problem, the Chandrasekhar 
algorithm was combined with F-reduction techniques [12] ,[17] to obtain 

reduced order approximations. 

A complete dynamic model for the system in which the elastic motions of 
the structure are coupled with the motions of the surrounding fluid results in 
a functional differential equation of neutral type (see Reference [6]). 
However, for this paper we shall use a simplified model based on the 
generalized Jones type approximations of the Wagner function described in 
[5]. The parameters used for the numerical examples below were obtained by 
applying the parameter identification scheme developed in [5] to experimental 
wind tunnel data. The resulting model is a five-dimensional delay- 

differential equation of the form 

x(t) = A^x(t) + AjX(t - r) + Bu(t) (15) 

x(s) = <p(s) ~r < s £ 0 (16) 


with output 


y(t) = Cx(t) 


(17) 
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In this model x = col(h,a,h,a,r) , where h is the plunge, a is the pitch 
angle and T represents a generalized aerodydnamic "lag state." The initial 
data was taken to be constant on [— r,0] . There is one control so that B is 
a 5 x 1 matrix and C = diag(c^ ,C 2 jC^ ,c^ ,c^) . The matrices Aq and 

are 5x5 with the only nonzero entry in A^ in the last row and last 
column. In particular, the time delay r = .05 and 


^ = 


-4.2106 

.5302 

1.0 

0.0 

-2.1187 


-31.2446 

-8.0098 

0.0 

1.0 

1262.16 


-4473.27 

563.315 

0.0 

0.0 

-2250.91 


-3704.37 

-7391.14 

0.0 

0.0 

-32863.5 


3.06111 

-1.15822 

0.0 

0.0 

-256.508 


A^^) = -47.00, AjCi.j) = 0, (i,j) * (5,5) 
B = [-81.6087 192.589 0.0 0.0 678.182] 1 
R = 10.0 

C = diag[ /50T0 /50T0 /IoTO /IoTO 1]. 
The initial state is the constant function 


<|>(s) = [-.80 .50 .055 .029 50.0] 
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and the control chosen represents a downward force applied at a point along 
the airfoil. All integrations were performed using a standard fourth order 
Runge-Kutta method with a fixed step size of h = .001. The Chandrasekhar 
algorithm was applied to the approximating system to obtain the gains K N (t). 
The resulting closed-loop system was integrated forward to calculate the 
optimal control and response. The system was initially solved on the 
interval [0, .25] . The results of the closed loop response (continuous line 
graph) and the unforced system (A graph) appear in Figures 2-6, and of 
particular note is that all closed loop responses approach zero as t + .25. 
In Reference [1] where a similar problem was treated, high frequency 
oscillations were obtained near t = .25 which may be due to numerical 
instabilities. 

The strength of the Chandrsekhar algorithm and F-reduction is revealed in 
a simple count of equations. The averaging scheme used divides [-r,0] 
into N subintervals and approximates the "history" of the equation by 

piecewise constant functions (see References [2] and [14] for details). 
Simulations of the above model were performed with N = 2,4,8,16,20 and 
convergence was obtained at N = 16. (The results in Figures 2-7 are for N = 
16). The averaging scheme results in an ordinary differential equation model 
that has 85 states. The gain, K^(t), is usually calculated by solving a 
matrix Riccati differential equation, and for N = 16 this necessitates 

solving 3655 equations. However, when the Chandrasekhar algorithm and F- 

reduction are applied, the number of states is reduced to 21, and it is 

necessary to solve only 126 equations to obtain K N (t). In this case, the 

calculation of the gain, state, and control took approximately 11 seconds on 
an IBM 3081. 
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We comment here that the Chandrasekhar equations were also integrated on 
the interval [0,1] to obtain the steady-state gain. In this case excellent 
performance was also obtained when the loop was closed. 




Figure 1 
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Figure 2 
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Figure 3 
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Figure 4 
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Figure 6 
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Figure 7 
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